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Abstract 

When the body is infected, it mounts an acute inflammatory response to rid 
itself of the pathogens and restore health. Uncontrolled acute inflammation due to 
infection is defined clinically as sepsis and can culminate in organ failure and death. 
We consider a three dimensional ordinary differential equation model of inflammation 
consisting of a pathogen, and two inflammatory mediators. The model reproduces the 
healthy outcome and diverse negative outcomes, depending on initial conditions and 
parameters. We analyze the various bifurcations between the different outcomes when 
key parameters are changed and suggest various therapeutic strategies. We suggest 
that the clinical condition of sepsis can arise from several distinct physiological states, 
each of which requires a different treatment approach. 
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1 Introduction 



The initial response of the body to an infection or trauma is called the acute inflam- 
matory response. This response is non-specific and is the first line of defense of the 
body against danger (?). It consists of a coordinated local and systemic mobilization 
of immune, endocrine and neurological mediators. In a healthy response, the inflam- 
matory response becomes activated, clears the pathogen (in the event of infection), 
begins a repair process and abates. However, inflammation itself can damage oth- 
erwise healthy cells which could then further stimulate inflammation. This runaway 
inflammation can lead to organ failure and death. Systemic inflammation accompa- 
nied by infection, based on its clinical manifestations, is defined as sepsis (?). Sepsis 
is a common and frequently fatal condition, with 750,000 cases annually in the United 
States alone in 1995 (?). 

Though much has been learned about the molecular and physiological pathways of 
the acute inflammatory response, this knowledge has not led to many effective ther- 
apies against sepsis. The sole approved drug therapy for severe sepsis is activated 
Protein-C, which only reduced mortality by 6% compared with controls in clinical 
trials (?; ?). One reason for the lack of effective treatments may be that the complex 
nature of the inflammatory response renders the effect of targeting isolated compo- 
nents of inflammation difficult to predict. Thus, mathematical modeling may provide 
insights into the global dynamics of the inflammatory process from which therapies 
may be developed. We propose that simple models of the acute inflammatory re- 
sponse can exhibit various outcomes and facilitate an understanding of the complex 
interactions between the various components of the response. 

We present a simple 3-dimensional model of the inflammatory response to infec- 
tion that captures the following clinically relevant scenarios: a healthy response where 
pathogen is cleared and the body returns to homeostasis, recurrent infection where 
inflammation is inadequate and the pathogen cannot be completely eliminated, per- 
sistent infectious inflammation where the pathogen levels and inflammation are high, 
persistent non-infectious inflammation where pathogen is cleared but inflammation 
persists and severe immuno- deficiency where pathogen has grown to saturation but 
the inflammatory response is very low. The model suggests that sepsis is a multi- 
faceted disease and narrowly targeted interventions are unlikely to succeed. 
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2 Reduced Model of Acute Inflammation 



Invading pathogens such as bacteria are rapidly detected by the body and an acute 
inflammatory response ensues. Among the first responders are phagocytic immune 
cells such as neutrophils and macrophages. These immune cells detect the bacterial 
cell components, become "activated" and begin to release pro-inflammatory cytokines 
such as Tumor Necrosis Factor-a (TNF-a), Interleukins (IL)-l, IL-6, IL-8 and High 
Motility Group Box-1 (HMGB-1) that activate more immune cells and recruit them 
to the sites of the infection (?). In addition, anti- inflammatory mediators such as IL- 
10 and Transforming Growth Factor-/? (TGF-/3) are also released which inhibit the 
production of the pro-inflammatory mediators. The activated phagocytic cells kill 
bacteria directly by engulfment and secretion of toxic chemicals such as oxygen free 
radicals (?). These substances can damage otherwise healthy cells. These damaged 
or dysfunctional cells can then induce more inflammation (?). 

Ideally, the inflammatory response eliminates the pathogen and then subsides. 
In some cases, the response might not be strong enough to clear the pathogen. In 
other cases, a positive feedback loop could arise between the early and late pro- 
inflammatory waves leading to a non-abating inflammatory response. Clinically, a 
sustained acute inflammatory response is manifested as septic shock and could cul- 
minate in organ failure and death. Though this is a simplification of the pathogenesis 
of sepsis, the idea is supported by the persistence of high levels of pro-inflammatory 
cytokines in non-survivors of sepsis (?; ?). 

A simplified picture of the acute inflammatory response is that an infectious 
pathogen triggers early pro-inflammatory responders which attempt to kill the pathogen. 
The early inflammatory mediators then activate later inflammatory mediators which 
can further excite the early mediators. This is the basis of our model which consists 
of three variables: 1) a pathogen p, which is an instigator of the innate immune 
response; 2) an early pro-inflammatory mediator m, which can be thought of as rep- 
resenting the combined effects of immune cells such as macrophages and neutrophils 
together with early pro-inflammatory mediators such as TNF-a and IL-1; 3) a late 
pro-inflammatory mediator I which represents a late pro-inflammatory feedback. This 
is a combined effect of cytokines such as IL-6, HMGB-1 and stimulatory effects of 
tissue damage and dysfunction. Although our model is extremely simple, it captures 
qualitatively all the salient features of a more complicated and biologically faithful 
model that is currently being developed (?). 

We consider mass-action type kinetics in a well-mixed volume. The dynamics 
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9 is an activation threshold and w is an activation width. All the variables and 
parameters are non-negative. The pathogen p obeys logistic growth and is killed 
when it interacts with m. The presence of p or / and m will stimulate the growth 
of m which also has an intrinsic death rate. This growth saturates as m increases 
towards unity mimicking the effects of cell depletion and anti-inflammatory cytokines. 
The late mediator / is recruited by m through a sigmoidal coupling function Q and is 
cleared with the rate of ki. In summary, the dynamics obey predator-prey dynamics 
with a delayed response. 

We note that the specific nature of the interactions of our model are not essential 
for the qualitative dynamics we find. As an example, the factor 1 — m in may 
be replaced by 1/(1 -|- m) or the form of the coupling function @ can be changed 
without affecting our conclusions. The important element is that there be an early 
and late inflammatory mediator with saturating kinetics. 

3 Dynamics of the Model 
3.1 Numerical examples 

The model © - ® exhibits behavior reminiscent of what is observed in clinical set- 
tings. Given an initial condition of p > 0, m > and / > 0, p grows, inducing m and / 
to grow and attempt to eliminate p. Depending on initial conditions and parameters, 
the ensuing orbits either approach stable fixed points or undergo oscillations, each 
having a physiological interpretation. We note from @ that a non-zero positive ini- 
tial value for m is necessary to generate an inflammatory response. The background 
level of late mediators is given by = {kim/ki)f{0) and is non-zero for our choice of 
parameters. Both these properties are consequences of our particular choice for the 
form of the model and not essential for the qualitative results we find. 

We must adopt one crucial element before we can interpret our results. In a 
healthy response to infection, the inflammatory response should become activated. 
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eliminate the bacteria and return to rest. However, in model (^3)-®) while p can 
diminish exponentially fast, it can never become completely eliminated in a finite 
time. Thus a consequence is that if m returns to rest after an inflammatory episode, 
p can re-grow from an infinitesimal quantity. This is an artifact of the model that 
arises because the approximation of p as a continuous variable breaks down when 
the population becomes small and the discreteness of the pathogen number becomes 
important. In this regime, a stochastic or agent based model where the pathogen can 
be completely eliminated is more appropriate. We finesse the discreteness problem 
by introducing a threshold for pathogen level. When p falls below this threshold, we 
consider it to be completely cleared. We propose that there is an effective thresh- 
old representing a single pathogen particle below which, on average, the pathogen 
population is eliminated. We will show how to calculate this threshold explicitly for 
various models in a future publication. 

We consider numerical examples for various values of kpm, kmp, and ki^- In 
Section 13.21 we show that these are natural bifurcation parameters of the system. 
The other parameters are fixed at kpg = 3, ki = 1, 9 = 1 and w = 0.5. We show the 
effects of varying 6 and w in Appendix IA.2.11 Numerical simulations and bifurcation 
plots were generated with XPPAUT (?). 

A healthy response to infection as seen in Fig. ^ corresponds to an orbit that 
spirals outwards so that p falls below threshold during the oscillation. The pathogen 
is then completely cleared which allows the inflammatory response to relax back to 
rest. 

In Fig. [2 the same parameters are used but the initial pathogen load is higher so 
that instead of returning to background levels, the inflammatory mediators are over- 
excited and remain elevated. We relate this situation to a state of persistent non- 
infectious inflammation where even though the pathogen is cleared, the inflammatory 
response does not abate. In our model, this state is a fixed point but in a real organism 
we expect that if this condition continued, it would eventually lead to multiple organ 
failure and death (?). 

Starting again from a healthy situation, if the pathogen susceptibility to the host's 
defenses {kpm) is decreased, we can enter a domain of persistent infectious inflam- 
mation where the inflammatory response is high but the pathogen still cannot be 
cleared as seen in Fig. |31 We would equate this condition with a severe septic state 
where both infection and inflammation are uncontrolled. In this case the damage 
caused by both the pathogen and inflammation are disrupting body function and if 
unabated death will result. Patients with systemic inflammation, with and without 
documented infection, are observed in clinical settings (?). 
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If we reduce the recruitment rate of / (kim), the healthy response can be turned 
into one of recurrent infection as seen in Fig. |^ (Changes in other parameters could 
also lead to recurrent infection as we see in the next section). In this case, low levels 
of infection persist indefinitely. This could be likened to infection with tuberculo- 
sis, yeast infections or low-grade bacterial infections that persist for long periods of 
time (?). Although, an organism could survive this state for a long duration it may 
eventually succumb. 

Finally, in Fig. |21 the activation rate of m due to p (kmp) is very weak, and we 
have a state of immunosuppression or immuno- deficiency where the pathogen grows 
to saturation and does not elicit any response from the body. This could happen if 
the immune-system had already been compromised by previous infection or trauma 
and then the organism is reinfected. Opportunistic bacterial and fungal infections 
have been observed in immuno-suppressed patient populations such as HIV infected 
patients, the elderly and those with organ transplants (?). These five scenarios are 
the only possible outcomes in the model. 



3.2 Fixed Points and Bifurcations 



These various regimes and transitions are best understood by examining the fixed 
points and associated bifurcations of the model Q-®- The fixed points satisfy the 
following conditions: 
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We consider three natural parameters a = kpm/kp, kmp, and b = hm/h, that appear 
in the fixed point conditions © and (jHJ. These three parameter combinations rep- 
resent the pathogen susceptibility to m compared to its growth rate (i.e. inverse of 
the pathogen virulence), the activation rate of early responders m due to p, and the 
effective recruitment rate of / due to m, respectively. 

The intersections of conditions © and (jSJ give the fixed points of the system. The 
dependence of the fixed points with the parameters is best observed as a projection 
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in the p — m plane as shown in Fig. Elwhere the fixed point conditions © and are 
plotted. Lines p = and m = are unaffected by changes in parameters. The line 



p = 1 — am 



(9) 



is affected only by a and sweeps across the p — m plane as a is decreased. The curve 



is affected only by the immune parameters b and kmp- It drops below p = as 6 is 
increased. Changing kmp affects the height and angle at which (fT(1|) intersects Q- 

There are five fixed points which we have labeled from FP 1 to FP 5. FP 1 is given 
by p = 0, m = 0, and / = {kim / h) f (0) . This fixed point is always unstable because 
the pathogen is growing. FP 2 is given by p = 1, m = 0, / = {kim/ki)f{0) and may 
be interpreted as a severely immuno- deficient state. Here the pathogen has grown to 
saturation, but there is no early immune response and the late response remains at 
the background level. This point is stable when the early immune response is very 
weak as shown in Appendix lA. 1.21 

FP 4 and FP 5 arise from a saddle node bifurcation when p = in H1U|) . This can 
be achieved by increasing b from zero as seen in the bifurcation plot in Fig. FP 4 
is always unstable and never represents any physiologically relevant scenario. When 
stable, FP 5 represents the persistent non-infectious inflammation fixed point. FP 3 
is given by the intersection of the line © with the curve (jlUj) . FP 3 could represent 
healthy, recurrent infection or persistent infectious inflammation states, depending on 
parameter values. Below, we vary the three parameter combinations a, kmp and b 
and examine the bifurcation plots. 

In the healthy scenario, a is large enough so that FP 3 is an unstable spiral and 
FP 5 is the only stable node {kpm = 20 in Fig.[7I). In this case, initial conditions close 
to FP 3 undergo oscillations that take the pathogen below the elimination threshold 
and are interpreted as healthy trajectories. A higher pathogen initial condition for 
the same set of parameters takes the trajectories directly to FP 5 and this behavior 
is interpreted as persistent non-infectious inflammation. When a is decreased, the 
slope of line © decreases and it sweeps through FP 4 rendering it unstable through 
a transcritical bifurcation. However, p is negative and thus FP 3 becomes unphysical 
{kpm = 10 in Fig.[7I). Hence FP 5 is the only stable fixed point for all initial conditions 
in this regime. As a is decreased further FP 3 crosses FP 5 and they exchange stability 
through another transcritical bifurcation. FP 3 then becomes the only stable fixed 
point and is the global attractor for all initial conditions. 
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As expected, starting from a healthy scenario, the outcome of the system becomes 
more and more severe with decreasing pathogen susceptibihty (i.e. increasing effec- 
tive pathogen virulence). If we were to increase susceptibility a starting from the 
healthy scenario, FP 3 changes into a stable spiral surrounded by an unstable limit 
cycle through a subcritical Hopf bifurcation. Now, trajectories within the unstable 
limit cycle spiral into the stable fixed point FP 3 resulting in recurrent infection and 
unstable spiral trajectories outside the unstable limit cycle will eventually dip below 
threshold and be interpreted as healthy. As FP 5 is still the attractor when the ini- 
tial p or m is very high, this set of parameters could lead to three different outcomes 
based on the initial conditions {kpm = 40 in Fig. [71). Increasing a further, increases 
the radius of the unstable limit cycle until it undergoes a homoclinic bifurcation when 
it collides with FP 1. Here FP 3 is a stable spiral and only recurrent infection is sup- 
ported. The model therefore predicts that pathogen clearance is not possible for a 
range of high pathogen susceptibility. However, as we increase susceptibility much 
further, FP 3 gets closer to the origin so that oscillations around it take p below 
threshold and may be interpreted as healthy. 

Varying k^p, the strength of the early response, does not affect the stability 
or position of FP 5, the persistent non-infectious inflammation point. Decreasing 
kmp takes a healthy state to a state where healthy and recurrent infection co-exist 
and finally to one of recurrent infection only as shown in Fig. |HJ Decreasing kmp 
below 1 — {kim/ki)f{0) makes FP 3 collide with FP 2 (the severely immuno-deficient 
state where p grows to saturation and m and / remain at background values). In a 
transcritical bifurcation, FP 2 becomes the stable fixed point of the system and FP 3 
becomes unphysical since m < 0. Increasing kmp from the healthy value reduces the 
value of p in FP 3 - this reduces the range of initial conditions for healthy behavior. 
However, for even very large values of kmp, FP 3 is an unstable spiral with complex 
eigenvalues supporting healthy behavior for initial conditions starting close to FP 3. 

Varying the strength of the late response, b affects both FP 3 and FP 5. As we 
decrease b, FP 3 undergoes a sub-critical Hopf followed by the homoclinic bifurcation 
of the unstable limit cycle similar to the above cases (Fig. |^. FP 4 and FP 5 
are created in a saddle-node bifurcation when ()1U() intersects p = by increasing 
b. However, there is also an upper limit to b beyond which healthy behavior is not 
supported. The curve (|1U() is below zero when it intersects @ and FP 3 is unphysical. 
Thus when the late response is too high, persistent non-infectious inflammation is the 
only possible outcome. 

In the preceding, although we have examined the behavior of the system by varying 
one parameter at a time, it gives us a picture of the global dynamics of the system. 
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Given the strong nonlinear saturation in the dynamics, orbits can either approach 
a fixed point, a hmit cycle, or a strange attractor. Given that FP 1 and FP 4 are 
always unstable and FP 2 is stable only when kmp is very low, this leaves FP 3 (or 
stable limit cycles around it) and FP 5 as the only candidates for global attractors. 
Ideally however, we would like FP 3 to be an unstable spiral or be surrounded by 
an unstable limit cycle, so that oscillations around it can be interpreted as healthy 
behavior. 

This gives an explanation for why healthy behavior always entails the risk of 
uncontrolled inflammation. In order for the pathogen to be cleared, we require an 
unstable spiral that can take p below threshold. However, if it were not for the (ex- 
ternally imposed) elimination threshold the orbit would eventually end up at some 
attractor and persistent non-infectious inflammation is the global attractor. Con- 
versely, if the curve (jlOj) has not yet intersected p = to form FP 4 and FP 5 (no 
persistent inflammation fixed point), FP 3 or a stable limit cycle around it (recurrent 
infection) is the global attractor. 

There is a possibility that FP 3 could undergo a supercritical Hopf bifurcation 
(discussed in Appendix lA. 1.3)) which leads to a stable limit cycle around it. How- 
ever, this is generically interpreted as recurrent infection because fine tuning would 
be required to ensure that the limit cycle takes p below threshold. Although, not 
possible in our model, a strange attractor might occur in higher dimensional models. 
We believe that this is unlikely given the large amounts of negative feedback and 
dissipation in the system. However, even if it were to exist, it may still not lead to 
a healthy situation as p may never be eliminated or the early and late inflammatory 
mediators could stay elevated albeit in a chaotic manner. 

Physiologically, this implies that in order to completely rid the body of a pathogen, 
the inflammatory response must respond strongly and remain elevated for a long 
enough time. If it responds too weakly, then the pathogen persists. If it abates too 
quickly, then there will be recurrent oscillatory infection. However, if it responds too 
strongly and too persistently then there is a risk that it will be self-sustaining even 
after the pathogen is cleared. Thus, there is a trade-off between being able to elimi- 
nate pathogens completely and risking nonabating inflammation. In Appendix lA. 21 
we calculate the parameter ranges where a healthy response is possible. 

4 Discussion 

Systemic inflammation and the ensuing organ damage is a major cause of mortality 
today (?; ?). This is also a disease created by modern medicine. Before the discovery 
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of advanced resuscitation techniques, patients could not be kept alive long enough 
for the condition to fully unfold. Patients often died from blood loss and severe in- 
fections before uncontrolled inflammation arose. With the advent and improvement 
of antibiotics and organ support therapy, the condition has become increasingly rel- 
evant (?). The incidence of systemic inflammation is also expected to increase with 
further advancement of medical technology and the aging of the population which 
are more susceptible (?). 

It is now recognized that the dysregulation of the underlying processes of acute 
inflammation can lead to multiple organ failure and death and that this is a com- 
mon pathway for diverse instigators such as trauma, hemorrhage and infection. Our 
simple model of acute inflammatory response to infection shows the various negative 
outcomes that arise from improper inflammatory response. These scenarios of persis- 
tent inflammation (with and without infection), immuno-suppression and recurrent 
infection have all been observed in critically ill patients (?; ?). Although strict cor- 
respondence with clinical reality is difficult to establish because of the simplicity of 
our model, we are able to classify our model behavior into the same broad categories. 
The model has a rich bifurcation structure and exploring it allows us to understand 
how changing parameters can take the system from one outcome to another. The 
model shows that in order to have a healthy response to infection the virulence of 
the pathogen cannot be too strong or too weak, the early pro-inflammatory response 
cannot be too weak, and the late response cannot be too strong or too weak. Analysis 
of this model gives clues to approach the problem of treating severe sepsis. Vastly 
different therapeutic strategies are called for to deal with the diverse negative out- 
comes. 

The strength of the acute inflammatory response varies in different individuals and 
may vary depending on age and environmental factors in the same individual. Studies 
suggest a strong genetic influence on the outcome of sepsis, and genetics may explain 
the wide variation in the individual response to infection (?). We examined the effect 
of the strength of early and late pro-inffammatory responses. The model suggests that 
only the strength of the late pro-inflammatory wave governs predisposition towards 
a state of persistent inflammation. No matter how exuberant the early wave may be, 
only controlling the feedback from the late wave can determine whether the outcome 
is healthy recovery or uncontrolled inflammation. Thus any therapy for persistent 
non-infectious inflammation must target the slow pro-inflammatory mediators. 

Experiments have found that survival was improved in infected mice when the 
moderately slow cytokine IL-6 was reduced (?). However, if too much IL-6 was 
removed then there was a detrimental effect. This result is consistent with our 
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model in that a small reduction lowered the possibility of a sustained inflamma- 
tory response but lowering it too much precluded the possibility of eliminating the 
bacteria. Evidence also has suggested that down-regulating HMGB-1, a late acting 
pro-inflammatory mediator might be a potential target for anti-sepsis therapies (?). 
Activated Protein- C which has been recently approved for treatment, is also partially 
a late anti- inflammatory agent (?). Previous therapeutic attempts have mostly fo- 
cused on down-regulating the early pro-inflammatory mediators and as predicted by 
the model have not shown great effectiveness (?; ?). 

On the other hand, if the patient is suffering from persistent infectious inflam- 
mation, then therapies must be aimed at both reducing the pathogen load and the 
late pro-inflammatory response. In this case, timing of the therapies may be impor- 
tant. It would be necessary to reduce the bacterial load first before reducing the 
inflammation. 

Conversely, low pathogen virulence or a weak immune response can lead to low 
level persistent or recurrent infection. Other theoretical models of infection have 
similar predictions. Persistence of the tuberculosis bacterium Mycobacterium tuber- 
culosis at low densities for extended periods has been suggested to be the result of 
slow growth rates (?). Similarly, a down-regulated immune response to Helicobacter 
pylori bacteria has been suggested to result in its persistent colonization of the hu- 
man stomach (?). In the clinical setting, patients in the ICU with decreased host 
defenses are susceptible to hospital- acquired infections (?). These infections which 
may be easily resolved in a healthy individual might result in unresolved inflammation 
and prove fatal to a compromised individual. Immuno-stimulatory therapy might be 
effective in such a situation. 

The model suggests that a healthy outcome is possible only when the risk of 
persistent inflammation is also present. A strong late immune response that increases 
the risk of unabated inflammation also ensures complete elimination of pathogen. 
This could explain why this risk has been retained by the evolution. 

In the clinical setting, inflammatory states are defined by symptoms and a few 
biological markers (?). For example, Systemic Inflammatory Response Syndrome 
(SIRS) is a condition characterized by - elevated respiratory and heart rates, fever, 
and an abnormal white blood cell count. The severity of sepsis is based on the 
presence of circulatory shock (low blood perfusion) and organ failure (?; ?). Similarly, 
treatment and support of these critically ill patients is largely based on clinical signs 
and a few biochemical and hematological parameters. However, the various outcomes 
of our model are based on the levels of the immune responders and the pathogen. 
Diverse physiological states as exemplifled in the model could be manifested by the 
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same clinical symptoms. Hence, one major problem for finding an effective treatment 
of sepsis is a diagnostic one. The model shows that sepsis should be considered as 
a set of distinct physiological disorders that require separate therapies even though 
they may have overlapping symptoms. Categorizing septic states based on levels of 
bio-markers rather than clinical symptoms would be the first step in addressing this 
problem and effort is ongoing in that direction (?; ?). 
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A Appendix 

A.l Eigenvalues of Fixed Points 

The Jacobian matrix of the model (^3) - @ is 
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A. 1.1 FP 1 
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This matrix is lower-triangular and the eigenvalues are kp, — 1, —ki. The fixed point 
is never stable as one of the eigenvalues kp is always positive. 



A. 1.2 FP 2 

For FP2: p = l,m = 0,1 = l^, the Jacobian is 
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The characteristic polynomial for the eigenvalue A is 

i-kp - X){-ki - \){krnp + /° - 1 - A) = 

with roots: X = —kp, X = —ki and A = k^p + 1^ — 1. Thus the condition for stability 
of FP2 is kmp < 1 — oi kmp < 1 — -^(1 — tanh(0/u;)). The other two eigenvalues 
are always negative. 



A. 1.3 FP 3 

FP3 is the solution of the following equations, 

p = 1 — am (14) 

kmpP ~l~ ^ 

/ = bf{m) (16) 
This may be reduced to a single transcendental equation for m: 

All the parameter combinations affect the position and eigenvalues of FP 3. Param- 
eter choices for which p or m of FP 3 is negative are unphysical. FP 3 could also 
have real or complex eigenvalues. FP 3 has complex eigenvalues and can support 
oscillatory behavior when ® intersects (\lU\i such that m of FP 3 is less than m of 
FP 4 (Fig. El). FP 3 can undergo Hopf bifurcations in this part of the phase-space 
- the subcritical Hopf bifurcations on varying the three parameter combinations are 
shown in Figs. [TllHl andlHl 

A supercritical Hopf bifurcation occurs when we change the coupling curve © so 
that 9 is small (low threshold for activation) and w is large (shallow coupling curve). 
In this case, for increasing b while keeping a and k^p fixed at moderate values, a state 
of persistent non-infectious inflammation will follow recurrent infection and a small 
zone of stable limit cycles. (Varying the other parameters, for small 9 and large w, 
does not result in Hopf bifurcations). In the case of small w (sharp coupling curve) 
and small 9 (low threshold of activation), it is possible that more than one fixed point 
satisfies the conditions for FP 3. In this case, curve H1U|1 intersects line © twice before 
dipping below zero. On varying the parameter combinations, recurrent infection is 
bistable with persistent inflammation and healthy behavior is never supported. These 
cases are discussed in greater detail in ADDendix lA.2.11 In our numerical examples, we 
have considered moderate w and large 9, which result in subcritical Hopf bifurcations. 
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A. 1.4 FP 4 and FP 5 



These fixed points are the solutions of the following equations, 



p = (18) 

I = ^f{m) (20) 

This may be reduced to one transcendental equation 

bfim) = (21) 
1 — m 

Equation (|21|) may have 0, 1 or 2 solutions depending on the parameters. Approxi- 
mate solutions can be found for H21() by substituting a piece-wise linear function such 
as /(m) = 0,m < 6 — w; f{m) = 1 + , 9 — w<m<6 + w; f{m) = 2,m > 6 + w. 
Using the above approximate function and assuming that FP 4 and 5 occur in region 
9 — w < m < 9 + w, we arrive at the following expressions for FP 4 and 5. 

p = (22) 

m = ^_ b_ ^23) 

I = bfim) (24) 

where q = 1 — 9 + w. 

This gives b > Aw/q^ as a lower limit on b for the existence of FP 4 and 5. When 
there are two solutions, they are formed through a saddle- node bifurcation as shown 
in Fig. FP 4 is a saddle point and FP 5 is a stable node. 

The eigenvalues of FP 4 and 5 can be found by substituting the solutions m = mP 
of (|21|) into the Jacobian (|11|). One of the eigenvalues is \i = kp — kpmin^ and the 
other two are functions of the parameters ki, kim, 9 and w. We note that k^p does 
not appear in the expressions for the position or eigenvalues of FP 4 and 5. 



A. 2 Regimes for Healthy Response 

Trajectories can be interpreted as healthy when FP 3 supports unstable oscillations 
as an unstable spiral or as a stable spiral surrounded by an unstable limit cycle. The 
unstable limit cycle is lost when it undergoes a non-generic homoclinic bifurcation 
when it collides with FP 1. The regimes of the various parameters (for fixed 9 and 
w) which maintain healthy behavior are discussed below. 

For healthy behavior, a should be large enough so that FP 3 is to the left of FP 5 
in Fig.Elor (1/a) < {q — \/ q^ — Aw/b)/2. Increasing kmp even over a large range, does 
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not make FP 3 unphysical, nor does FP 3 lose its complex eigenvalues - thus even 
for very large k^p, healthy oscillations are supported. The upper limit for a and the 
lower limit for kmp to maintain healthy oscillations are difficult to find in closed form 
as these variations result in the radius of the limit cycle increasing and oscillations 
being lost through a non-generic homoclinic bifurcation. 

The strength of the late immune response b should not be too high or too low. If b 
is too high, then FP 3 has p < and hence healthy behavior cannot be supported. In 
order for the late immune response to not be too high, we require that p of condition 
(|1U|) be positive so that FP 3 is not unphysical. Applying this to (|1U|) gives the 
condition 

^ < l-tanh(0/u;)- ^^^^ 
If b is too low, FP 4 and 5 are not yet formed and (jlOj) is well above p = so FP 3 
remains the global attractor and is a stable spiral. To avoid too low of a response, we 
would like FP 4 and 5 to exist so that unstable oscillations from FP 3 are eventually 
attracted to FP 5. Thus the condition for existence for FP 4 and 5, as calculated 
above in lA.1.4l is 

4:W 

A. 2.1 Effect of Varying the Coupling Curve 

The shape of the coupling curve (0]) (i.e. parameters 6 and w) can alter the results. 
We used 9 = 1 and w = 0.5 in our analysis. Parameter 9 sets the threshold where / 
is activated by m and parameter w gives the steepness of recruitment. 

From numerical simulations, we find that when w is too small, FP 3 never under- 
goes a Hopf bifurcation, and as ki^/ki is increased, recurrent infection and persistent 
inflammation are bistable (For example, w = 0.005 with 'healthy' parameter set in 
Fig. ^) . As discussed in Appendix IA.1.31 small w could result in a shallow enough 
curve (|10|) that intersects line @ more than once. The two fixed points replacing 
FP 3 arise in a saddle-node bifurcation. The stable node is bistable with FP 5 imply- 
ing that recurrent infection is bistable with persistent non-infectious inflammation. 
When w is too large {w = 1), FP 3 is a stable spiral that becomes unphysical as 
kim/ki increases. Thus the system supports recurrent infection followed by persistent 
inflammation and healthy behavior does not exist. In summary, for w too small or 
large the subcritical Hopf bifurcation of FP 3 (and hence healthy behavior) does not 
occur for any of the parameter combinations. 

When 9 is small (9 = 0.1), FP 3 undergoes a supercritical Hopf bifurcation as b is 
varied but this cannot be generically interpreted as healthy behavior as oscillations 
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about a stable limit cycle may not always take p below threshold. As 9 is increased, 
subcritical Hopf bifurcations are possible. The range of parameter b, given by H25|) 
and H26|) . where oscillatory behavior is observed also increases as 9 is increased. Thus, 
large 9 and moderate w would maximize the region for healthy response. 

In the above analysis, we have looked at varying parameter combinations rather 
than individual parameters for conciseness of the discussion. Varying the parameters 
such as kp and kpm individually does not change the bifurcation diagrams qualitatively 
and our interpretations of the effect of pathogen susceptibility remain unchanged. 
Similarly, varying and ki individually does not change our conclusions about the 
strength of late response. The phase-space analyzed is restricted between < m < 1 
since condition (|1U() prevents m from increasing beyond one. Using any other function 
such as 1/(1 + m) to saturate m can change the upper bound of m but the same 
behavior is retained. 
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Figure 1: Time courses of p (black), m (red) and / (blue) for the healthy response. We show 
orbits with the p elimination threshold (solid line) and without (dotted line). With the 
threshold (set at po = 0.0005), when p drops below pq, p is set to zero and m and / return 
to background values. Without the threshold, the orbits spiral outwards. Parameters used 
are kpm = 30, kp = 3, kmp = 25, kim = 15, ki = 1. Initial conditions are p(0) = 0.01, m(0) = 
0.05, /(O) = 0.539. 



17 




Time 



Figure 2: In persistent non-infectious inflammation, p (black) is eliminated but m (red) and 
/ (blue) remain elevated. Parameters used are kpm = 30, kp = 3, kmp = 25, kim = 15, ki = 1. 
Initial conditions are p{0) = 0.2, m(0) = 0.05, /(O) = 0.539. 
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Figure 3: In persistent infectious inflammation, p (black) cannot be eliminated and m 
(red) and / (blue) remain elevated (as in severe sepsis). Parameters used are kpm = 3,kp = 
3, kmp = 25, kim = 15, ki = 1. Initial conditions are p{0) = 0.01, m(0) = 0.05, /(O) = 0.539. 
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Figure 4: In recurrent infection, p (black) is low and m (red) and / (blue) remain above 
background values but not very high (as in a low-grade infection). Parameters used are 
kpm = 30, fcp = 3,kmp = 25, kim = 5,ki = 1. Initial conditions are p{0) = 0.01, m(0) = 
0.05, /(O) = 0.179. 
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Figure 5: In severe immuno- deficiency^ p (black) grows to saturation and m (red) and / 
(blue) remains low or absent. Parameters are kpm = 30, kp = 3, kmp = 0.4, kim = 15, ki = 
l,p(0) = 0.01, m(0) = 0.05, /(O) = 0.539. 



20 



0.2 0.4 0.6 0.8 1 

M (Early Mediators) 

Figure 6: Projection of fixed-point conditions onto the p-m plane. Projected directions of 
field-fiow are also marked. The red curves correspond to conditions p = and and 
the black curves correspond to conditions m = and (fTUI) . The fixed points are at the 
intersections of black and red curves and the direction of fiow along a curve reverses at 
fixed points. The projection of the phase-space for a typical set of parameters is shown. 
FP 1 and FP 4 are saddles. FP 3 is a spiral which is stable or unstable based on the actual 
values of the parameters. FP 5 is a stable fixed point. FP 2 (not shown) is also a saddle 
in this case. 
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Figure 7: Bifurcation plot showing outcome as kp^ and hence a = kp^/kp is increased. For 
low kpm, FP 3 interpreted as persistent infectious inflammation is the only fixed point of the 
system. As kpm is increased to ~ 2.7, FP 5 {persistent infectious inflammation) becomes 
the attractor of the system. FP 3 becomes unphysical and moves leftward in the p-m plane 
(see inset). At kpm ~ 11, FP 3 undergoes a transcritical bifurcation with FP 4, becoming 
physical again. It is an unstable spiral, oscillations around which are interpreted as healthy. 
On further increasing kpm, FP 3 undergoes a subcritical Hopf bifurcation and is surrounded 
by an unstable limit-cycle, oscillations around which are interpreted as healthy. Trajectories 
within the limit-cycle, which spiral into FP 3, are interpreted as recurrent infection. Healthy 
behavior is lost to a homoclinic bifurcation at kpm ~ 57. For large kpm, FP 3 and FP 5 
are bistable (i.e. recurrent infection and persistent non-infectious inflammation are possible 
outcomes of the system). (Inset : In the p-m plane, line sweeps across the plane causing 
the above bifurcations.) 
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Figure 8: Bifurcation plot showing outcome increased. For very low kmp{< 0.5), 

FP 2 is stable and is interpreted as the severely immuno- deficient state. As k^p is increased, 
FP 3 becomes stable through a transcritical bifurcation and recurrent infection is a possible 
outcome of the system. As k^p is increased further FP 3 undergoes a subcritical Hopf 
bifurcation {k^p ~ 22). Oscillations around the unstable limit cycle surrounding FP 3 are 
interpreted as healthy and are lost through a homoclinic bifurcation {kmp ~ 9). Trajectories 
within the limit-cycle spiral into FP 3 and are interpreted as recurrent infection. Beyond 
the Hopf bifurcation {k^p > 22), FP 3 remains an unstable spiral, up to very large kmp 
{kmp ~ 2000) and oscillations around FP 3 are interpreted as healthy. On varying kmp, 
only the position and eigenvalues of FP 3 are changed. FP 5 remains unaffected which 
is interpreted as persistent infectious inflammation and remains a possible outcome of the 
system. If curve fllOj) is shallow enough (6 was small enough) so that FP 5 is not yet created, 
then FP 3 remains the stable global attractor through changes in kmp- (Inset : In the p-m 
plane, FP 3 varies in position and stability as (fTUI) varies, causing the above bifurcations. 
FP 5 remains unaffected.) 
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Figure 9: Bifurcation plot showing outcome as kim or b = kim/h is increased. For low fc/^, 
FP 3 is a stable spiral representing recurrent infection. As kim is increased, FP 4 and 5 
are created through a saddle-node bifurcation at kim ~ 7 and the persistent non-infectious 
inflammation is now a possible outcome. At kim ~ 14, FP 3 undergoes a subcritical Hopf 
bifurcation and is surrounded by an unstable limit cycle. The oscillations around the 
unstable limit cycle are interpreted as healthy behavior and trajectories inside the limit 
cycle spiral into FP 3 and are interpreted as recurrent infection. Healthy behavior is lost at 
kim ~ 10 due to a homoclinic bifurcation. Beyond the Hopf bifurcation, healthy oscillations 
are still supported by FP 3 which is now an unstable spiral. At kim ~ 21, FP 3 meets 
FP 4 and becomes unphysical through a transcritical bifurcation (as its p < 0). Beyond 
this point, FP 5 is the only stable attractor and persistent non-infectious inflammation is 
the only possible outcome of the system. (Inset : In p-m plane, (fTUI) descends as kim is 
increased causing the above bifurcations.) 
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